Two serendipitous low-mass LMC clusters discovered with HST. 
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ABSTRACT 

We present V and I photometry of two open clusters in the LMC down to V ^ 26. The 
clusters were imaged with the Wide Field and Planetary Camera - ^ on board of the Hubble 
Space Telescope, as part of the Medium Deep Survey Key-Project. Both are low luminosity 



lO '• {My ~ —3.5), low mass systems (M ~ IO^Mq). The chance discovery of these two clusters 
in two parallel WFPC2 fields suggests a significant incompleteness in the LMC cluster census 

,—1 ■ near the bar. One of the clusters is roughly elliptical and compact, with a steep light profile, a 

^ ' central surface brightness nvi^) ~ 20.2 mag/arcsec^, half-light radius rhi ~ 0.9 pc (total visual 

. major diameter D ~ 3 pc) and an estimated mass M ~ 15OOM0. ^From the colour-magnitude 

'"p_^ diagram and isochrone fits we estimate its age as r ~ 2 — 5 x 10^ years. Its mass function has a 

^ , fitted slope of F = Alog(f){M) / AlogM = -1.8 ±0.7 in the range probed (0.9 ~ M/Mq ~ 4.5). 

^ The other cluster is more irregular and sparser, having shallower density and surface brightness 

^ _ profiles. We obtain F = — 1.2±0.4, and estimate its mass as M ~ 4OOM0. A derived upper 

^ " limit for its age is r ~ 5 x 10^ years. Both clusters have mass functions with slopes similar 



to that of R136, a massive LMC cluster, for which HST results indicate F ~ —1.2. They also 
seem to be relaxed in their cores and well contained in their tidal radii. 
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1 INTRODUCTION 



The Large Magellanic Cloud (LMC) contains a vast number of star clusters, with ages 
varying from 10^ to 10-^° years. The LMC cluster system is thus suitable for studying the 
evolution of physical properties of clusters such as mass, radius, light and density profiles, 
luminosity functions (LFs) and mass functions (MFs). These latter may depend not only on 
age, but also on metal abundance or environment and provide relevant information to the 
physics of star formation from fragmenting gas clouds (McClure et al. 1986, Larson 1991, 
1992). 

Derivation of relevant physical parameters in LMC clusters, as well as in other dense stel- 
lar systems, has been limited by lack of spatial resolution. This situation improved considerably 
with the refurbishment of the Hubble Space Telescope (HST), allowing crowding problems in 
dense stellar systems to be substantially reduced (Elson et al. 1995, de Marchi & Paresce 
1995a, b, Hinitcr ct al. 1996, Santiago et al. 1996). However, the faint end of the cluster lu- 
minosity function in the LMC has not yet been targeted by deep photometric studies. HST 
imaging of clusters and associations published so far has concentrated on Galaxy clusters or on 
rich systems in the Local Group. Almost nothing is known about small LMC clusters, whose 
detection is often difficult and whose properties are harder to determine observationally. The 
luminosity function of LMC clusters seems to rise steeply in the low luminosity domain (Elson 
& Fall 1985). In fact, the total number of detected LMC clusters has been steadily increasing 
and is currently believed to be about ~ 4500 (Hodge 1988). 

In this paper, we report observations of two low luminosity {My ~ —3.5) clusters in the 
LMC. They were detected in two parallel Wide Field and Planetary Camera - 2 (WFPC2) 
images as part of the Medium Deep Survey HST key-project (MDS). This chance discovery 
of two clusters suggests that such objects may be very common in the LMC. If our two MDS 
fields are typical, the implied surface density would be considerably larger than that inferred 
by Elson & Fall (1985) or by the deep photographic study of Hodge (1988), at least in the 
vicinities of the LMC bar or the 30 Doradus region. 

In §2, we present HST photometry in two bands for the two MDS fields; we show colour- 
magnitude diagrams (CMDs) and discuss completeness corrections and photometric calibration 
issues. We then use the data to extract surface brightness and stellar density profiles for the 
clusters, and to estimate their ages (§3). In §4, we compute their LFs and MFs. We briefly 
discuss their dynamical state as well. In §5, we discuss the results and present our conclusions. 

2 THE DATA 

2.1 Data Reduction 

The two LMC fields studied in this work were imaged in parallel mode with WFPC2, as 
part of the Medium Deep Survey (MDS) key project. Field 1 is located at a = 05 : 35 : 36.8 
and 5 = -69 24 23 (J2000). Field 2 is located at a = 05 36 52.3 and 5 = -69 37 59.6 (J2000). 
Both fields are at the eastern end of the LMC bar, southwest of 30 Doradus, close to NGC2050 
and NGC2048, respectively. Each field contains one of the small clusters reported in this work. 
We hereafter refer to the cluster located in field 1 (field 2) as Cl (C2). 
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Field 1 was imaged with two 500s exposures using the HST F814W (I) filter and two 500s 
exposures with the HST F606W (V) filter. For field 2, two exposures in each filter were taken 
as weU: 4900s and 1100s in F814W and 4900s and 1000s in F606W. 

The raw data were processed with the standard pipeline procedure, which corrects for 
instrumental effects (Holtzman et al. 1995a). The two exposures in each field/filter combination 
were then coadded and median filtered; the lower instrumental value was used at each pixel 
position to eliminate cosmic rays. Field 2 F814W frames were offset by about 3" from each 
other and were registered to a common position before coadding. 

2.2 Sample Selection 

An object list was obtained separately for each field/filter configuration from the final 
coadded image. The IRAF DAOPHOT package (Stetson 1987) was used for this as well as for 
aperture and psf fitting photometry. We adopted a detection threshold of 3 cr, where a is the 
standard deviation in the background counts of each chip. We worked only with the 3 Wide 
Field Camera - 2 (WFC2) chips, since the Planetary Camera (PC) chip would not significantly 
increase the sample size. 

Inspection of the images showed that most of the objects detected were real stars. How- 
ever, some spurious detections occurred, especially around bright stars. In order to clean up 
the sample, a point spread function (psf) template was created from a few bright unsaturated 
stars in each WFC2 chip and then fitted to all remaining objects using ALLSTAR. This task 
gives a a-nd a sharpness parameter (s). The dependence of these parameters on Veoe mag- 
nitude for chip 4, field 1 is shown in Figure 1. There are two loci in the sharpness diagram: 
inspection of the images revealed that the clump of objects with s > 0.1 contains predomi- 
nantly spurious detections (also attested by their large x^), whereas the low sharpness locus 
was almost entirely made up of stars. An additional cut in < 2 was also applied to eliminate 
the few remaining objects whose fit to the psf template was not satisfactory. Also, in order 
to avoid saturation effects, all objects brighter than /814 = 18 or Vgoe = 19 were eliminated. 
Figure 1 is typical of the other two WFC2 chips in field 1, both for F606W and F814W. 

Field 2 has a longer exposure time, leading to enhanced crowding and more saturated 
stars. The saturation magnitudes measured for field 2 are Igid = 20.5 and Vqqq = 21.5. 
Registration of the F814W images reduced the quality of the psf fits. The s x /814 and 
X -^814 diagrams did not display as clear boundaries between stellar and non-stellar objects. 
Thus, star selection for field 2 was restricted to the F606W image, although object detection 
was carried out independently in both filters. The faint saturation limits in field 2 made it 
difficult to select bright stars based on the x ^606 and s x Veoe diagrams. Thus, only 
stars with Veoe ~ 21.5 were selected this way. Brighter stars were selected from a smoothed 
version of the coadded images, obtained by applying a 3 pixel (~ 0.3") Gaussian. This filter is 
wide enough to eliminate hot pixels and psf features, but narrow enough to allow the detection 
of bright objects. 

The number of stars was typically 6000-7000 per chip in field 1 and 8000 in field 2. Of 
these, about 5000 had both IgiA and Veoe magnitudes available. 

2.3 Photometry 
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Even though sample selection was mostly based on psf fitting, the magnitudes and colors 
used in the analysis came from aperture photometry inside a radius r = 2 pix (r = 0.2"). Vqqq 
and /814 magnitudes measured this way led to a narrower CMD than that based on the psf 
magnitudes. The choice of radius is a compromise between the need to bypass centering and 
undersampling problems and the need to avoid light contamination from neighbouring sources. 
An aperture correction of 0.24 mag had to be applied to both HST filters in order to account 
for the light outside the aperture (Holtzman et al. 1995a). 

For field 2, saturation prevented magnitudes and colors of upper main sequence and the 
red giant branch (RGB) stars to be measured. In order to bypass this problem, aperture 
photometry for the bright stars in field 2 (Veoe < 21.5) was carried out in the F814W and 
F606W frames with the shortest exposure time. These have much less stringent saturation 
limits (JgM ~ 18.75, Vgoe ~ 19.75). 

Finally, an attempt was made to measure magnitudes for field 2 stars with /814 < 18.75 or 
Vqoq < 19.75 using only the pixels within the circular ring between 2 and 3 pixel radius. Again 
the shortest exposures were used. Given the shape of the WFC2's psf, this in principle allows 
one to push saturation levels towards brighter magnitudes. The quality of the magnitudes 
obtained with this procedure was tested with several isolated, bright and unsaturated stars. 
Only approximate magnitudes could be measured, the uncertainty being ~ 0.35 mag. About 
50 such bright stars were added to the sample in each field 2 chip. 

2.4 Photometric Calibration 

The data were calibrated to the Johnson- Cousins system using the "synthetic" transfor- 
mation equations listed in Table 10 of Holtzman et al. (1995b). 

Reddening corrections were determined by comparing isochrones and the aperture cor- 
rected instrumental CMDs until the observational and theoretical main-sequences and RGBs 
coincided. We used Yale isochrones for that purpose (Green et al. 1987), and assumed a dis- 
tance modulus of m — M = 18.5 to the LMC (Panagia et al. 1991). The Yale isochrones had 
first to be converted to the HST photometric system by applying the inverse of the calibration 
equations listed by Holtzmann et al. (1995b). For field 1, a -^(Veoe — -^814) = 0.12 was obtained. 
For field 2 E{Vqoq — /814) = 0.17 provided the best fit. Given the difficulty in fitting all the 
CMD features, however, these values may be uncertain by as much as 0.1 mag. 

The reddening and aperture corrected HST colour-magnitude diagrams (CMDs) for all 
stars in fields 1 and 2 are shown in Figure 2. Also shown are Yale isochrones for Z=0.01 stars 
with ages of 0, 200, 500, 1000 and 2000 x 10^ years (Myrs). Only non-saturated stars are 
included. The CMD for field 1 (panel 2a) shows a main sequence ranging from Veoe ~ 18.5 
down to Vqoq ~ 26. The RGB is also clearly visible in both panels, but especially for field 1, 
where the shorter exposure time and the more accurate sample selection and photometry allow 
even the red giant clump (Veoe ^ 18-5 — 19) to stand out. The larger number of objects in 
Field 2 which lie outside the main-sequence and RGB reflects the limitation in star selection, 
which for field 2 was based on the F606W image only. 

The isochrones nicely fit the main sequences of both panels justifying the adopted red- 
dening and distance modulus. In particular, the isochrones bracket the width of the upper 
main sequence, which should be made up of fairly young stars and for which the choice of 
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metallicity is justified (Olszewski et al. 1991). A detailed analysis and discussion of the field 
stellar populations and star formation history in these fields was left to Elson et al. (1997). 

2.5 Completeness Corrections 

Completeness functions were obtained independently for field and cluster stars, since they 
are known to depend on crowding. Completeness functions were measured mostly for WFPC2 
chip 4, field 1 and WFPC2 chip 3, field 2, where CI and C2 are respectively located. 

For the field stars, we ran a total of 40 reahzations of the DAOPHOT.ADDSTAR task 
for each HST filter, 5 realizations for each of 8 magnitude bins, spanning the range 20 < 
/814, Vqoq < 28. In each realization, 200 stars were added to randomly chosen sections of the 
original frame, each section being 400x400 pixels wide and situated away from the cluster. The 
section images containing the artificial stars were then put through the same detection and 
star classification processes as the real data. Only artificial stars whose input and observed 
magnitudes were within 0.3 mag of each other were considered as detected. The average 
fraction, out of the 5 realizations for each magnitude bin, of artificial objects detected and 
classified as stars was taken as the completeness value at that magnitude. 

For the cluster regions we made 40 realizations with 10 artificial stars for each of 12 mag- 
nitude bins within the range 20 < Isi^, Vqqq < 26. Experiments showed that CI completeness 
does not depend strongly on position within the cluster region. This is consistent with the 
compact and steeply declining profile for this cluster (§3). For C2, two completeness functions 
were assigned, one for its core region and the other for the outskirts. 

The Veoe completeness functions are shown in Figure 3. Cluster completeness falls more 
rapidly with magnitude than field completeness due to more severe crowding. The differences 
in completeness between fields 1 and 2 are almost always smaller than the error bars. This 
applies both to cluster and field stars completeness functions. 

A faint cut-off limit was applied to the data in order to avoid large shot-noise errors. 
Magnitude limits of /814 = 25 (/814 = 24) and Vqqq = 26 {Vqqq = 25) were applied to field 
(cluster) stars. The magnitude errors (1 a) at the cut-off limits are SIsia ~ 0.2 for /814 = 25 and 
SVqog ~ 0.3 for Veoe — 26. These limits were used in the derivation of structural parameters, 
density and surface brightness profiles, LFs and MFs presented in §3 and §4. 

The /814 completeness functions behave similarly to the ones shown in Figure 3. Since the 
final sample used in this work is that made up of objects with both a Vqqq and a /814 magnitude, 
a joint completeness function has computed. The additional incompleteness caused by the 
requirement of V (I) band detection was quantified by simply multiplying, at each magnitude 
level, the I (V) completeness function, by the fraction of sources detected (and classified as 
stars) in I (V) which made into the final sample. 

3 THE LOW-MASS LMC CLUSTERS: MORPHOLOGY, STRUCTURE 
AND AGE 

3.1 Morphology and Structure 

Vqqq band images of CI and C2 are shown in Figures 4 and 5. CI is located at a = 
05 35 35.1 and 6 ^ -69 23 48.8, whereas C2's centre is at a = 05 36 51.9 and 6 - -69 38 12.1 
(J2000). These positions are about 20' and 30' away from the center of the 30 Dor region, 
respectively. 
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CI has visually determined major and minor diameters of 12" x 10". Assuming a 
distance modulus of m — M = 18.5 for the LMC, this corresponds to D x d ~ 3 x 2.5 pc 
(1 pc = 4.1"). 

In Figure 6 we show the stellar number density (panel a) and the surface brightness 
profiles (panel b) for CI. The upper curve on both panels gives the profile uncorrected for 
contamination by field stars. The error bars include Poisson fluctuations as well as uncertainties 
in the completeness functions. The dots show background corrected profiles and incorporate 
the additional statistical uncertainty associated with the subtracted background stars. Arrows 
indicate upper limits. In this case the tip of the arrow is at the most probable value and its 
upper end corresponds to the la deviation from this value. The horizontal lines in both panels 
indicate the background levels, determined separately for each WFC2 chip (dotted lines), and 
linearly interpolated into the cluster region (solid line). This latter was used as the best 
estimate of the field contamination at the cluster position. The profiles shown include only 
stars in the range —0.5 < My < 6.5. 

A clear excess of stars is visible out to i? ~ 2 pc (~ 80 pixels), beyond which the stellar 
number density merges with that estimated for the LMC field (panel 6a) . There are at least as 
many cluster stars as background ones within this radius. Cl's central regions have a roughly 
constant surface brightness, nv ~ 20.2 mag/arcscc^ (panel 6b). Beyond R ~ 0.5 pc, however, 
the surface brightness profile falls steeply with radius, flattening out again at -R ~ 2 pc. This 
outer extension is not visually noticeable (see flg. 4) and is not present in the number density 
proflle either (panel 6a). It could be an artifact caused by underestimation of background 
surface brightness levels; as panel 6b itself shows, the background fleld pv varies by some 0.3 
mag with position within fleld 1. However, the uncorrected nv proflle is still brighter than the 
highest background level estimated from the WFC2 chips. Since CI lies in an intersection of 
several associations and star forming regions, the excess of surface brightness beyond i? ~ 2 pc 
may be caused by these larger scale structures in the LMC. We return to this issue in §4.1. 

The stellar density and surface brightness proflles for C2 are shown in Figure 7. Again 
only stars with —0.5 < My < 6.5 contribute to the proflles. Because C2 is closer to the chip 
border than Cl, its proflles do not extend as far from the cluster center as in the case of CI. 
C2 has visual diameters of 19" x 16" {D x d = 4.6 x 3.9 pc). Its visual appearance suggests 
a sparser and more irregular cluster showing some substructure. In spite of the low stellar 
number density contrast relative to the background, panel 7a shows a systematic excess of 
stars out to R ^ 1.5 pc, although this excess is hardly signiflcant beyond ~ 1 pc. C2's light 
proflle (panel 7b), on the other hand, is above that of the contaminating fleld out to ~ 1.5 pc, 
the excess brightness relative to the background being still signiflcant all the way out to the 
edge of the chip. We get //y(0) ~ 20.3 mag/arcsec^ for C2, comparable to Cl. Its /ly proflle 
is shallower than that of Cl. 

The inferred structural parameters for Cl and C2 are listed in Table 1, including sizes 
and central densities. 

3.2 Ages 

The top panels of Figure 8 show the CMD for stars within boxes of 16" (4 pc) on a side 
centered on Cl (panel a) and on C2 (panel b). A total of 268 and 310 stars are shown in 
panels a and b, respectively. The data are in the Johnson- Cousins system and are corrected 
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for aperture and extinction effects, as described in §2.3 and §2.4. The dashed hnes indicate 
saturation levels. They correspond to a fixed Vqqq cut-off (see §2.2 and §2.3). The lower panels 
show CMDs for stars in the neighbourhood of each cluster, for comparison. The entire area 
outside the cluster in the WFC2 chip where it is located was used as comparison field. The field 
CMDs shown include a randomly selected fraction of the field stars, so that the numbers of 
cluster and background stars in any region of the CMD can be directly compared. They clearly 
differ in the upper main sequence: there are 16 main sequence stars with V < 20 m panel a, 8 of 
which have V < 19. In the corresponding comparison field (panel c), these numbers are 8 and 
1. Similarly, only 3 main sequence stars have F < 20 in panel d, whereas in the corresponding 
cluster CMD there are 14 such stars. Thus, most of the upper main sequence stars are real 
cluster members and we expect field contamination not to afi^ect age estimates from isochrone 
fitting. 

Yale isochrones corresponding to Z=0.01 stars with ages of 0, 200, 500, 1000 and 2000 
Myrs are shown in the upper panels. This chosen metallicity is typical of LMC clusters (Ol- 
szewski et al. 1991). Its associated uncertainty {5[Fe/H] ~ 0.15) has a smaller effect on the 
isochrone fits than the reddening and saturation effects. The 500 Myrs isochrone is the one 
that best fits Cl's upper main sequence. However, a few saturated stars exist within CI and 
have been left out of Figure 8. The presence of stars brighter than the 500 Myrs turn-off would 
indicate a younger age. In fact, adjustments in the amount of extinction or in the metallicity 
assumed for the stars would allow an age as low as r = 200 Myrs for CI. 

Assigning an age to C2 is harder, given the larger photometric errors and stronger sat- 
uration effects. Another problem is that the observed main sequence is a bit redder than the 
isochrones, suggesting a larger reddening within C2's region than elsewhere in field 2. Adjust- 
ing the reddening values in order to match the theoretical and observational MSs and taking 
into account the presence of several saturated stars within C2, we can only set an upper limit 
of T ~ 500 Myrs to the age of C2. 

The derived ages for CI and C2 are also quoted in Table 1. In the next two sections we 
determine the luminosity and mass functions for both clusters and evaluate their slopes. 

4 Luminosity and Mass Functions 

4-1 Cluster Luminosity Functions 

In Figure 9 we show completeness corrected luminosity functions for CI and C2. Field 
contamination was eliminated by subtracting the field LF from that within the cluster region. 
The field LF was taken to be the average over several control regions equidistant from the 
cluster. The cluster regions used for determining their LFs were circles of 8" radius centered 
on each cluster. 

CI seems to have a slightly steeper LF than C2. LF slopes were obtained from linear fits 
to the points. We obtained 7 = Alog^{Mv)/AMv = 0.19 ± 0.03 for CI in the range My < 6 
and 7 = 0.12 ± 0.05 for C2 in the range My < 4. The best fit lines are shown in the figure. 
Upper limits were not included in the fits. Given the small number of LF bins and the low 
contrast of the clusters, the slope differences are not significant. In fact, the two LF slopes 
are similar in the common range used for the fits {My < 3). C2 LF basically ends beyond 
that while CI LF steepens. The 7 values are in agreement with that inferred from the work of 
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Flower et al. (1980) for NGC1868, a rich LMC cluster with similar age CI and C2 but larger 
mass. On the other hand, the inferred values for 7 are smaller than those typically fitted to 
younger globular clusters and to stellar associations in the LMC (Vallenari et al. 1993, Will et 
al. 1995a,b), 

In Figure 10, we show CI LFs for 3 radial bins; the two first are 1 pc wide, the last is 2 
pc wide. All 3 LFs shown are field subtracted and were scaled to the entire cluster area. The 
number of stars per unit area decreases with radius, as expected. The outermost ring contains 
only 3 bins in My with numbers significantly above the background. This ring includes the 
stars that make up the excess surface brightness seen beyond R — 2 pc in panel 6b. Despite 
the uncertainties, the LF looks shallower in panel 10c than in 10b, consistently with figure 6, 
where only an excess of light, not stars, is seen. In the inner rings, on the other hand, the LF 
becomes steeper with radius, providing evidence for mass segregation within CI. Thus, the 
bright stars beyond J? ~ 2 pc are likely to be background stars, belonging either to the general 
LMC field or to some stellar association superposed to CI. In fact, CI is situated in a rather 
messy border region between different clusters and associations, among them NGC2050, LH96, 
DEM261 and NGC157 (see catalogs by Lucke & Hodge 1970, Davies et al. 1976). The best fit 
slopes for the LFs in the two inner rings are 7 = 0.09 ± 0.04 and 7 = 0.25 ± 0.02. The outer 
ring lacks enough points for reliable fits to be made. 

It was not possible to split C2's LF into radial sectors, given its smaller contrast with the 
background. 

4.2 Cluster Mass Functions 

Comparing the MFs of clusters with different masses, ages, metallicities or environments 
can contribute to the debate about the universality of the initial mass function (IMF) and its 
evolution. 

The MF can be derived from the LF with the aid of a mass-luminosity (M-L) relation, 
which usually depends on metallicity and is often uncertain. However, the M-L relation is 
reasonably well known for the typical metallicity of the LMC, Z ~ 0.008 {[Fe/H] ~ —0.3), in 
the range of luminosities covered in this work. The mass functions for CI and C2 are shown 
in Figure 11. The M-L conversion was based on the Yale isochrone that best fits the clusters 
CMDs. Field subtraction proceeded in the same way as with the LFs. For field stars we used 
the same M-L relation as for the clusters. The derived mass of the evolved stars will obviously 
be in error, but these are subtracted off with the field, having little or no effect on the cluster 
MFs shown. 

Linear fits were made to the data points leading to F = Alog(f){M)/ AlogM = —1.8 ± 0.7 
and —1.2 ± 0.4 for Cl and C2, respectively. For C2, however, a single power-law fit to the MF 

is inappropriate, since its MF seems to be steeper for M ~ 2 Mq. 

The derived slopes are subject to several sources of error. Isochrones with different 
metallicities or based on stellar models which incorporate convective core overshooting would 
change the M-L relation and therefore the MF slopes. The effect, however, is known to be 
small: ~ 0.1 in F (Elson et al. 1989, Sagar et al. 1991). Unresolved binaries lead to an observed 
MF which is fiatter than the single star one. Sagar & Richtler (1991) investigated this issue 
and concluded that the amplitude of the effect is a function of the binary fraction and the 
MF slope itself. The observed slopes derived in this paper may be consistent with a true 
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r ~ — 2 and —1.5 for CI and C2, respectively, if the fraction of binaries is as large as 0.5. 
Finally, completeness errors may distort the actual shapes of LFs and MFs, especially if mass 
segregation is present, as seems to be the case in CI. Civen the small size of the clusters, it was 
impossible to split them into many annuli for which completeness and mass functions could be 
measured independently. Given all these uncertainties we cannot rule out a single MF shape 
accounting for both of them. 

The LF and MF slopes are listed in Table 2, along with estimates of total luminosities 
{My), masses (M) and mass-to-light ratios (M/L). These latter were obtained by adding up 

the contributions of all non-saturated stars with My ~ 6.5 (M ~ 0.9 Mq). Thus, the absolute 
magnitude and mass estimates listed should be considered as lower limits. 

The resulting MF slopes are bracketed by most found in the literature. These were all 
derived for clusters much larger than the ones studied here. Elson et al. (1989), for instance, 

obtained —1.8 ~ F ~ — 0.8 for six rich young (r ~ 30 — 50 Myrs) LMC clusters in the mass 

range 1.5 ~ M/Mq ~ 6. Such flat IMFs were in large contrast with the results of Mateo 

(1988), who found slopes (—2.5 ~ F ~ — 4.6) for six LMC clusters spanning a wider range of 
ages and metallicities. Sagar & Richtler (1991) derived intermediary slopes, F ~ —2.1, for five 
LMC clusters within the mass range 2 < M/Mq < 14 and spanning an age range of 10 — 100 
Myrs. 

More recently and using HST data. Hunter et al. (1995) obtained F = —1.22 for R136, 
a ~ 2 X 10^ Mq, very young and compact system in the center of 30 Doradus, in the range 
2.8 < M/Mq < 15. Our results suggest flat MF slopes for low mass LMC clusters as well. 

A more meaningful comparison would be with open clusters in the Galaxy. Ground- 
based studies include that of Sagar et al. (1986) who derived MF slopes for 11 open clusters. 
One of them (N1778C) has an age comparable to those of CI and C2 in §3.2. Its slope 
is F = —2.04 ± 0.44 in a mass range comparable to ours. This is consistent with CI and 
marginally consistent with C2. Other more recent ground-based works have led to shallower 
slopes, F ~ — 1, but usually in the mass range M ~ 0.5 Mq (Comeron et al. 1996, Hambly et 
al. 1995 and references therein). 

4-3 Internal dynamics 

We now use the integrated properties and structural parameters listed in Tables 1 and 
2 to assess the importance of dynamical effects such as two body relaxation and tidal forces. 
Knowledge of the internal dynamical state of CI and C2 may help determining the extent to 
which such low mass clusters could contribute to the field population. 

Using the LMC tidal field derived by Elson, Fall & Freeman (1987), we estimate CI and 

C2 to have similar tidal radii, in the range 6 ~ (pc) ~ 15, implying rt > r^is (see Table 1). 
Therefore, the clusters are so far not tidally truncated. Both, however, may have experienced 
stronger tidal fields or had close encounters with other more massive clusters or associations 
during their ~ 10^ years of existence. 

Is two-body relaxation relevant to the internal dynamics of CI and C2? The relaxation 
time scale for stars with a mass M at a distance r from the centre is given by 



2 X 10^ 



(2 pc) ( 



0.2 kms-i 



V 




(1) 
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where v is the typical velocity of such a star within the cluster. Assuming hydrostatic 
equilibrium in the central regions, the central dispersion velocity in the cluster is given by 

<j\Q)r.Gporl (2) 

where Tc is the core radius and po is the central mass density. Using the relation between 
projected and spatial densities given by Djorgovski (1993), we estimate po ~ 16 Mq/pc^ for 
CI and po ~ 6 Mq/pc^ for C2. These values for the central densities imply (t(0) ~ 0.4 km 
for CI and cr(0) ~ 0.2 km s~^ for C2. Inserting these velocities into equation (1) we infer that 

trei ~ 100 Myrs for a 1 Mq star in the central regions of both clusters. Therefore, both CI 
and C2 should be relaxed in their central parts. The evidence for mass segregation in CI is 
consistent with that. Dispersion velocities just about twice the values inferred would disrupt 
the clusters in ~ 50 Myrs. Given their estimated ages, CI and C2 are probably bound. 

5 DISCUSSION AND CONCLUSIONS 

We presented V and I photometry of two WFPC2 fields located near the eastern end of 
the LMC bar. Stars as faint as My ~ 7 (M ~ O.SM©) were detected in each of them. Each 
field contained one small open cluster, implying a large number density of such systems in the 
LMC. 

CI is roughly symmetrical in shape, has a steep density profile and a mass function slope 
ofr = — 1.8±0.7. ^From isochrone fits, we infer an age of r ~ 200 — 500 Myrs for it. Its 
estimated mass is ~ 1500 M© and its absolute magnitude. My ~ —3.5. The derived values for 
the luminosity, mass, MF slope and age are mutually consistent in light of the recent stellar 
population synthesis models of Girardi et al. (1995). These authors use the photometric models 
of single stellar populations calculated by Bertelli et al. (1994) in order to revise the relation 
between integrated photometric properties, age and metallicity of LMC clusters. ^From their 
Figure 13, we infer that a M ~ 2 x lO^M© cluster with a MF slope a bit shallower than a 
Salpeter one (F = —2.35; Salpeter 1955) would have My ~ —3.5 at an age r ~ 10^ years. 

Marginal evidence of mass segregation within r ^ 2 pc from Cl's centre was found. This 
is consistent with the short relaxation time scale expected for M ~ 1 Mq stars in its central 

regions. Tidal effects from the LMC should not be relevant; we derive rhi ~ 1 from Cl's 
surface brightness profile, which places its member stars well inside the estimated tidal radius 
(8 ~ rt (pc) ^ 15). Given the central density and size estimates for Cl, its stars should have a 
maximum central velocity dispersion of ~ 0.4 km s~^. It would be interesting to confirm that 
with observations. 

C2 seems less massive but just as luminous as Cl (M ~ 4OOM0, My ~ —3.5). This 

implies fiatter luminosity and mass functions. We obtain F = —1.2 ±0.4 for C2. For My ~ 3.5, 
C2's LF is similar to Cl's but it drops off for fainter magnitudes. C2 is more irregular than 
Cl and has shallower density and surface brightness profiles. We could only derive an upper 
limit of T ~ 500 Myrs to the age of C2 based on its CMD due to saturation and reddening 
uncertainties. Based on the results of Girardi et al. (1995) and our mass. My and F estimates, 

we would obtain more stringent limits: r ~ lOOMyrs. C2's size was harder to quantify, given 
its fiat profile. However, just as Cl, it is also likely to be contained within its tidal radius and 
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to have undergone significant core relaxation: trei ~ 100 Myrs for M ~ IMq stars within the 
central 2 pc. 

The chance discovery of two small clusters in two MDS fields within the LMC suggests 
that these could be more common in the LMC bar or 30 Doradus region than previously 
anticipated. Such systems may have significantly contributed to the field star population if 
they were even more common in the past and got disrupted by larger clusters or by the LMC 
tidal field. Even though that does not seem to be the case for CI and C2, other similar systems 
may have been subjected to stronger tidal fields from both the LMC or its bar or from more 
massive clusters or associations. Alternatively many such clusters might have been unbound 
since their birth, their member stars streaming away after a few hundred Myrs. It would 
interesting to confirm or not the existence of a large population of low mass and luminosity 
clusters in the LMC with other observations using the high resolution of HST. 
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Table 1. Clusters structural parameters. 



Parameter 


CI 


C2 


(Jo (stars/pc^) 


30 ±8 


8±6 


Ra (pc) 


~ 2 


- 1.5 


Hvi^) (mag/arcsec^) 


20.2 ±0.4 


20.3 ±0.5 


Rf^ (pc) 


~ 2 


~ 2 


Rhi (pc) 


~ 0.9 


~ 0.7 


Dyis ^ dyis (pc) 


2.9 X 2.4 


4.6 X 3.9 


T{Myrs) 


200-500 


~ 500 
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Table 2. LF and MF slopes. 



Parameter CI C2 



7 = Alog^Mv)/AlogMv 0.19 ±0.03 0.12 ±0.05 

F = Mog^{M)/MogM -1.8 ±0.7 -1.2 ±0.4 

3.5 3.5 

Mass^ (Mo) - 1500 - 400 

M/L (Mq/Lq) ~ 0.7 ~ 0.2 



1 Due to stars in the range -0.5 ^ My ^ 6.5 (0.9 ~ M/Mq ~ 4.5) only. 
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FIGURE CAPTIONS 

1- The sharpness (s) and parameters, as obtained by fitting each detected object in chip 
4, field 1 to the WFC2 psf, are shown as a function of Vqqq magnitudes, a) s. b) x^- 

2- a) CMD for all 3 WFC2 chips in field 1. The magnitudes and colours are in the HST 
WFPC2 system. The data are corrected for aperture and reddening effects. Superposed 
are Yale isochrones corresponding to Z=0.01 stars with ages 0, 200, 500, 1000 and 2000 
Myrs. b) The same as in a but now for the 3 WFC2 chips in field 2. 

3- a) Vqoq completeness functions for field stars in field 1, chip 4 (squares and dashed line) 
and in field 2, chip 3 (crosses and dotted line). Error bars are standard deviations from 
different realizations of ADDSTAR, as explained in the text, b) Same as in a but now 
cluster completeness functions are shown: squares and dashed line: CI; crosses and dotted 
lines: inner and outer C2. Error bars are shown only for one of the cluster completeness 
functions to avoid confusion. 

4- WFPC-2 image of LMC field 1 cluster, CI. The whole field shown is about 50" on a side. 

5- WFPC-2 image of LMC field 2 cluster, C2. The whole field shown is about 40" on a side. 

6- a) Stellar density profile within the region occupied by CI. The upper curve gives the 
surface density of stars as a function of distance from the cluster centre, without any 
correction for field contamination. The points indicate the background corrected profile. 
Arrows correspond to upper limits. The adopted background density level is given by the 
solid horizontal line. The dotted horizontal lines show the background densities estimated 
separately for each WFC2 chip. 6^ V band surface brightness profiles and background 
levels corresponding to the number density ones shown in panel a. 

7- a) The same as in panel 6a but now for C2. b) The same as in panel 6b but now for C2. 

8- a) The colour magnitude diagram of stars within CI region. Yale isochrones of 0, 200, 
500, 1000 and 2000 Myrs have been superposed to the data. The dashed line indicates the 
saturation level, b) The same as in panel a but now for C2. c) CMD for field comparison 
region close to CI. d) CMD for field comparison region close to C2. 

9- The luminosity functions of CI and C2 are shown as solid and open squares, respectively. 
The error bars incorporate Poisson fiuctuations as well as uncertainties in the complete- 
ness corrections and in the background subtraction. Upper limits are shown with arrows. 
The lines show power-law fits to the data. 

10- The luminosity function for CI in different radial annuli, as indicated. The numbers 
in all panels are scaled to the total cluster area. Error bars again incorporate statistical 
uncertainties as well as completeness and background errors. Arrows correspond to upper 
limits. The solid line shown in the two upper panels is the best fit to a power-law. 

11- Mass functions of CI (solid symbols) and C2 (open symbols). The error bars incorporate 
Poisson fluctuations as well as uncertainties in the completeness corrections and in the 
background subtraction. The lines show the best fit power-laws in each case (CI: solid 
line; C2: dashed line). 
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ABSTRACT 

We present V and I photometry of two open clusters in the LMC down ioV^ 26. The 
clusters were imaged with the Wide Field and Planetary Camera - 2 on board of the Hubble 
Space Telescope, as part of the Medium Deep Survey Key-Project. Both are low luminosity 
{My ~ —3.5), low mass systems (M ~ IO^Mq). The chance discovery of these two clusters 
in two parallel WFPC2 fields suggests a significant incompleteness in the LMC cluster census 
near the bar. One of the clusters is roughly elliptical and compact, with a steep light profile, a 
central surface brightness //^/(O) ^ 20.2 mag/arcsec^, half-light radius Vhi ^ 0.9 pc (total visual 
major diameter D ~ 3 pc) and an estimated mass M ~ 1500Mq. From the colour-magnitude 
diagram and isochrone fits we estimate its age as r ~ 2 — 5 x 10^ years. Its mass function has a 
fitted slope of P = Alog(l){M) / AlogM = -1.8 ±0.7 in the range probed (0.9 ^ M/Mq ^ 4.5). 
The other cluster is more irregular and sparser, having shallower density and surface brightness 
profiles. We obtain P = —1.2 ± 0.4, and estimate its mass as M ~ 4OOM0. A derived upper 
limit for its age is r 5 x 10^ years. Both clusters have mass functions with slopes similar 
to that of R136, a massive LMC cluster, for which HST results indicate P ~ —1.2. They also 
seem to be relaxed in their cores and well contained in their tidal radii. 
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1 INTRODUCTION 



The Large Magellanic Cloud (LMC) contains a vast number of star clusters, with ages 
varying from 10^ to 10^° years. The LMC cluster system is thus suitable for studying the 
evolution of physical properties of clusters such as mass, radius, light and density profiles, 
luminosity functions (LFs) and mass functions (MFs). These latter may depend not only on 
age, but also on metal abundance or environment and provide relevant information to the 
physics of star formation from fragmenting gas clouds (McClure et al. 1986, Larson 1991, 
1992). 

Derivation of relevant physical parameters in LMC clusters, as well as in other dense stel- 
lar systems, has been limited by lack of spatial resolution. This situation improved considerably 
with the refurbishment of the Hubble Space Telescope (HST), allowing crowding problems in 
dense stellar systems to be substantially reduced (Elson et al. 1995, de Marchi & Paresce 
1995a, b. Hunter et al. 1996, Santiago et al. 1996). However, the faint end of the cluster lu- 
minosity function in the LMC has not yet been targeted by deep photometric studies. HST 
imaging of clusters and associations published so far has concentrated on Galaxy clusters or on 
rich systems in the Local Group. Almost nothing is known about small LMC clusters, whose 
detection is often difficult and whose properties are harder to determine observationally. The 
luminosity function of LMC clusters seems to rise steeply in the low luminosity domain (Elson 
& Fall 1985). In fact, the total number of detected LMC clusters has been steadily increasing 
and is currently believed to be about ~ 4500 (Hodge 1988). 

In this paper, we report observations of two low luminosity (My ^ —3.5) clusters in the 
LMC. They were detected in two parallel Wide Field and Planetary Camera - 2 (WFPC2) 
images as part of the Medium Deep Survey HST key-project (MDS). This chance discovery 
of two clusters suggests that such objects may be very common in the LMC. If our two MDS 
fields are typical, the implied surface density would be considerably larger than that inferred 
by Elson & Fall (1985) or by the deep photographic study of Hodge (1988), at least in the 
vicinities of the LMC bar or the 30 Doradus region. 

In §2, we present HST photometry in two bands for the two MDS fields; we show colour- 
magnitude diagrams (CMDs) and discuss completeness corrections and photometric calibration 
issues. We then use the data to extract surface brightness and stellar density profiles for the 
clusters, and to estimate their ages (§3). In §4, we compute their LFs and MFs. We briefiy 
discuss their dynamical state as well. In §5, we discuss the results and present our conclusions. 

2 THE DATA 

2.1 Data Reduction 

The two LMC fields studied in this work were imaged in parallel mode with WFPC2, as 
part of the Medium Deep Survey (MDS) key project. Field 1 is located at a = 05 : 35 : 36.8 
and 6 = -69 24 23 (J2000). Field 2 is located at a = 05 36 52.3 and 5 = -69 37 59.6 (J2000). 
Both fields are at the eastern end of the LMC bar, southwest of 30 Doradus, close to NGC2050 
and NGC2048, respectively. Each field contains one of the small clusters reported in this work. 
We hereafter refer to the cluster located in field 1 (field 2) as CI (C2). 
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Field 1 was imaged with two 500s exposures using the HST F814W (I) filter and two 500s 
exposures with the HST F606W (V) filter. For field 2, two exposiu^es in each filter were taken 
as well: 4900s and 1100s in F814W and 4900s and 1000s in F606W. 

The raw data were processed with the standard pipeline procedure, which corrects for 
instrumental effects (Holtzman et al. 1995a). The two exposures in each field/filter combination 
were then coadded and median filtered; the lower instrumental value was used at each pixel 
position to eliminate cosmic rays. Field 2 F814W frames were offset by about 3" from each 
other and were registered to a common position before coadding. 

2.2 Sample Selection 

An object list was obtained separately for each field/filter configuration from the final 
coadded image. The IRAF DAOPHOT package (Stetson 1987) was used for this as well as for 
aperture and psf fitting photometry. We adopted a detection threshold of 3 a, where a is the 
standard deviation in the background counts of each chip. We worked only with the 3 Wide 
Field Camera - 2 (WFC2) chips, since the Planetary Camera (PC) chip would not significantly 
increase the sample size. 

Inspection of the images showed that most of the objects detected were real stars. How- 
ever, some spurious detections occurred, especially around bright stars. In order to clean up 
the sample, a point spread function (psf) template was created from a few bright unsaturated 
stars in each WFC2 chip and then fitted to all remaining objects using ALLSTAR. This task 
gives a a sharpness parameter (s). The dependence of these parameters on Veoe mag- 

nitude for chip 4, field 1 is shown in Figure 1. There are two loci in the sharpness diagram: 
inspection of the images revealed that the clump of objects with s > 0.1 contains predomi- 
nantly spurious detections (also attested by their large x^), whereas the low sharpness locus 
was almost entirely made up of stars. An additional cut in < 2 was also applied to eliminate 
the few remaining objects whose fit to the psf template was not satisfactory. Also, in order 
to avoid saturation effects, all objects brighter than J814 = 18 or Veoe = 19 were eliminated. 
Figure 1 is typical of the other two WFC2 chips in field 1, both for F606W and F814W. 

Field 2 has a longer exposure time, leading to enhanced crowding and more saturated 
stars. The saturation magnitudes measured for field 2 are /814 = 20.5 and Feoe = 21.5. 
Registration of the F814W images reduced the quality of the psf fits. The s x /8i4 and 
X /814 diagrams did not display as clear boundaries between stellar and non-stellar objects. 
Thus, star selection for field 2 was restricted to the F606W image, although object detection 
was carried out independently in both filters. The faint saturation limits in field 2 made it 
difficult to select bright stars based on the ^ Veoe and s x Veoe diagrams. Thus, only 
stars with Veoe ^ 21.5 were selected this way. Brighter stars were selected from a smoothed 
version of the coadded images, obtained by applying a 3 pixel (~ 0.3") Gaussian. This filter is 
wide enough to eliminate hot pixels and psf features, but narrow enough to allow the detection 
of bright objects. 

The number of stars was typically 6000-7000 per chip in field 1 and 8000 in field 2. Of 
these, about 5000 had both Isi4 and Veoe magnitudes available. 

2.3 Photometry 
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Even though sample selection was mostly based on psf fitting, the magnitudes and colors 
used in the analysis came from aperture photometry inside a radius r = 2 pix (r = 0.2"). Veoe 
and /814 magnitudes measured this way led to a narrower CMD than that based on the psf 
magnitudes. The choice of radius is a compromise between the need to bypass centering and 
undersampling problems and the need to avoid light contamination from neighbouring sources. 
An aperture correction of 0.24 mag had to be applied to both HST filters in order to account 
for the light outside the aperture (Holtzman et al. 1995a). 

For field 2, saturation prevented magnitudes and colors of upper main sequence and the 
red giant branch (RGB) stars to be measured. In order to bypass this problem, aperture 
photometry for the bright stars in field 2 (Feoe < 21.5) was carried out in the F814W and 
F606W frames with the shortest exposure time. These have much less stringent saturation 
limits (/814 ^ 18.75, Veoe ^ 19.75). 

Finally, an attempt was made to measure magnitudes for field 2 stars with /814 < 18.75 or 
1^606 < 19.75 using only the pixels within the circular ring between 2 and 3 pixel radius. Again 
the shortest exposures were used. Given the shape of the WFC2's psf, this in principle allows 
one to push saturation levels towards brighter magnitudes. The quality of the magnitudes 
obtained with this procedure was tested with several isolated, bright and unsaturated stars. 
Only approximate magnitudes could be measured, the uncertainty being ^0.35 mag. About 
50 such bright stars were added to the sample in each field 2 chip. 

2.4 Photometric Calibration 

The data were calibrated to the Johnson-Cousins system using the "synthetic" transfor- 
mation equations listed in Table 10 of Holtzman et al. (1995b). 

Reddening corrections were determined by comparing isochrones and the aperture cor- 
rected instrumental CMDs until the observational and theoretical main-sequences and RGBs 
coincided. We used Yale isochrones for that purpose (Green et al. 1987), and assumed a dis- 
tance modulus ofm — M=18.5to the LMC (Panagia et al. 1991). The Yale isochrones had 
first to be converted to the HST photometric system by applying the inverse of the calibration 
equations listed by Holtzmann et al. (1995b). For field 1, a -^(Veoe — -^814) = 0.12 was obtained. 
For field 2 ^'(Veoe ~ -^814) = 0.17 provided the best fit. Given the difficulty in fitting all the 
CMD features, however, these values may be uncertain by as much as 0.1 mag. 

The reddening and aperture corrected HST colour-magnitude diagrams (CMDs) for all 
stars in fields 1 and 2 are shown in Figure 2. Also shown are Yale isochrones for Z=0.01 stars 
with ages of 0, 200, 500, 1000 and 2000 x 10^ years (Myrs). Only non-saturated stars are 
included. The CMD for field 1 (panel 2a) shows a main sequence ranging from Veoe ^ 18.5 
down to Veoe ^ 26. The RGB is also clearly visible in both panels, but especially for field 1, 
where the shorter exposure time and the more accurate sample selection and photometry allow 
even the red giant clump (Veoe ^ 18.5 — 19) to stand out. The larger number of objects in 
Field 2 which lie outside the main-sequence and RGB refiects the limitation in star selection, 
which for field 2 was based on the F606W image only. 

The isochrones nicely fit the main sequences of both panels justifying the adopted red- 
dening and distance modulus. In particular, the isochrones bracket the width of the upper 
main sequence, which should be made up of fairly young stars and for which the choice of 
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metallicity is justified (Olszewski et al. 1991). A detailed analysis and discussion of the field 
stellar populations and star formation history in these fields was left to Elson et al. (1997). 

2.5 Completeness Corrections 

Completeness functions were obtained independently for field and cluster stars, since they 
are known to depend on crowding. Completeness functions were measured mostly for WFPC2 
chip 4, field 1 and WFPC2 chip 3, field 2, where CI and C2 are respectively located. 

For the field stars, we ran a total of 40 reahzations of the DAOPHOT.ADDSTAR task 
for each HST filter, 5 realizations for each of 8 magnitude bins, spanning the range 20 < 
-^814) 1^606 < 28. In each realization, 200 stars were added to randomly chosen sections of the 
original frame, each section being 400x400 pixels wide and situated away from the cluster. The 
section images containing the artificial stars were then put through the same detection and 
star classification processes as the real data. Only artificial stars whose input and observed 
magnitudes were within 0.3 mag of each other were considered as detected. The average 
fraction, out of the 5 realizations for each magnitude bin, of artificial objects detected and 
classified as stars was taken as the completeness value at that magnitude. 

For the cluster regions we made 40 realizations with 10 artificial stars for each of 12 mag- 
nitude bins within the range 20 < /8i4) 1^606 < 26. Experiments showed that CI completeness 
does not depend strongly on position within the cluster region. This is consistent with the 
compact and steeply declining profile for this cluster (§3). For C2, two completeness functions 
were assigned, one for its core region and the other for the outskirts. 

The Veoe completeness functions are shown in Figure 3. Cluster completeness falls more 
rapidly with magnitude than field completeness due to more severe crowding. The differences 
in completeness between fields 1 and 2 are almost always smaller than the error bars. This 
applies both to cluster and field stars completeness functions. 

A faint cut-off limit was applied to the data in order to avoid large shot-noise errors. 
Magnitude limits of /814 = 25 (/814 = 24) and Veoe = 26 (Veoe = 25) were applied to field 
(cluster) stars. The magnitude errors (1 a) at the cut-off limits are 5/314 ^ 0.2 for /8i4 = 25 and 
5^606 ^ 0.3 for Veoe = 26. These limits were used in the derivation of structural parameters, 
density and surface brightness profiles, LFs and MFs presented in §3 and §4. 

The /814 completeness functions behave similarly to the ones shown in Figure 3. Since the 
final sample used in this work is that made up of objects with both a Veoe and a /814 magnitude, 
a joint completeness function has computed. The additional incompleteness caused by the 
requirement of V (I) band detection was quantified by simply multiplying, at each magnitude 
level, the I (V) completeness function, by the fraction of sources detected (and classified as 
stars) in I (V) which made into the final sample. 

3 THE LOW-MASS LMC CLUSTERS: MORPHOLOGY, STRUCTURE 
AND AGE 

3.1 Morphology and Structure 

Veoe band images of CI and C2 are shown in Figures 4 and 5. CI is located at a = 
05 35 35.1 and S = -69 23 48.8, whereas C2's centre is at a = 05 36 51.9 and S = -69 38 12.1 
(J2000). These positions are about 20' and 30' away from the center of the 30 Dor region, 
respectively. 
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CI has visually determined major and minor diameters of 12" x 10". Assuming a 
distance modulus of m — M = 18.5 for the LMC, this corresponds toD x d ~ 3 x 2.5 pc 
(1 pc = 4.1"). 

In Figure 6 we show the stellar number density (panel a) and the surface brightness 
profiles (panel b) for CI. The upper curve on both panels gives the profile uncorrected for 
contamination by field stars. The error bars include Poisson fluctuations as well as uncertainties 
in the completeness functions. The dots show background corrected profiles and incorporate 
the additional statistical uncertainty associated with the subtracted background stars. Arrows 
indicate upper limits. In this case the tip of the arrow is at the most probable value and its 
upper end corresponds to the la deviation from this value. The horizontal lines in both panels 
indicate the background levels, determined separately for each WFC2 chip (dotted lines), and 
linearly interpolated into the cluster region (solid line). This latter was used as the best 
estimate of the field contamination at the cluster position. The profiles shown include only 
stars in the range —0.5 < My < 6.5. 

A clear excess of stars is visible out to i? ~ 2 pc (~ 80 pixels), beyond which the stellar 
number density merges with that estimated for the LMC field (panel 6a). There are at least as 
many cluster stars as background ones within this radius. Cl's central regions have a roughly 
constant surface brightness, jiy ^ 20.2 mag/arcsec^ (panel 6b). Beyond R ~ 0.5 pc, however, 
the surface brightness profile falls steeply with radius, fiattening out again at ^ 2 pc. This 
outer extension is not visually noticeable (see fig. 4) and is not present in the number density 
profile either (panel 6a). It could be an artifact caused by underestimation of background 
surface brightness levels; as panel 6b itself shows, the background field fiy varies by some 0.3 
mag with position within field 1. However, the uncorrected jiy profile is still brighter than the 
highest background level estimated from the WFC2 chips. Since CI lies in an intersection of 
several associations and star forming regions, the excess of surface brightness beyond R ^ 2 pc 
may be caused by these larger scale structures in the LMC. We return to this issue in §4.1. 

The stellar density and surface brightness profiles for C2 are shown in Figure 7. Again 
only stars with —0.5 < My < 6.5 contribute to the profiles. Because C2 is closer to the chip 
border than CI, its profiles do not extend as far from the cluster center as in the case of CI. 
C2 has visual diameters of 19" x 16" {D x d = 4.6 x 3.9 pc). Its visual appearance suggests 
a sparser and more irregular cluster showing some substructure. In spite of the low stellar 
number density contrast relative to the background, panel 7a shows a systematic excess of 
stars out to ^ 1.5 pc, although this excess is hardly significant beyond i? ~ 1 pc. C2's light 
profile (panel 7b), on the other hand, is above that of the contaminating field out to R ^ 1.5 pc, 
the excess brightness relative to the background being still significant all the way out to the 
edge of the chip. We get fiv{0) ^ 20.3 mag/arcsec^ for C2, comparable to CI. Its /iy profile 
is shallower than that of CI. 

The inferred structural parameters for CI and C2 are listed in Table 1, including sizes 
and central densities. 

3.2 Ages 

The top panels of Figiu-e 8 show the CMD for stars within boxes of 16" (4 pc) on a side 
centered on CI (panel a) and on C2 (panel b). A total of 268 and 310 stars are shown in 
panels a and b, respectively. The data are in the Johnson-Cousins system and are corrected 
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for aperture and extinction effects, as described in §2.3 and §2.4. The dashed lines indicate 
saturation levels. They correspond to a fixed Veoe cut-off (see §2.2 and §2.3). The lower panels 
show CMDs for stars in the neighbourhood of each cluster, for comparison. The entire area 
outside the cluster in the WFC2 chip where it is located was used as comparison field. The field 
CMDs shown include a randomly selected fraction of the field stars, so that the numbers of 
cluster and background stars in any region of the CMD can be directly compared. They clearly 
differ in the upper main sequence: there are 16 main sequence stars with < 20 in panel a, 8 of 
which have V < 19. In the corresponding comparison field (panel c), these numbers are 8 and 
1. Similarly, only 3 main sequence stars have V" < 20 in panel d, whereas in the corresponding 
cluster CMD there are 14 such stars. Thus, most of the upper main sequence stars are real 
cluster members and we expect field contamination not to affect age estimates from isochrone 
fitting. 

Yale isochrones corresponding to Z=0.01 stars with ages of 0, 200, 500, 1000 and 2000 
Myrs are shown in the upper panels. This chosen metallicity is typical of LMC clusters (Ol- 
szewski et al. 1991). Its associated uncertainty {5[Fe/H] ~ 0.15) has a smaller effect on the 
isochrone fits than the reddening and saturation effects. The 500 Myrs isochrone is the one 
that best fits Cl's upper main sequence. However, a few saturated stars exist within CI and 
have been left out of Figure 8. The presence of stars brighter than the 500 Myrs turn-off would 
indicate a younger age. In fact, adjustments in the amount of extinction or in the metallicity 
assumed for the stars would allow an age as low as r = 200 Myrs for CI. 

Assigning an age to C2 is harder, given the larger photometric errors and stronger sat- 
uration effects. Another problem is that the observed main sequence is a bit redder than the 
isochrones, suggesting a larger reddening within C2's region than elsewhere in field 2. Adjust- 
ing the reddening values in order to match the theoretical and observational MSs and taking 
into account the presence of several saturated stars within C2, we can only set an upper limit 
of r ^ 500 Myrs to the age of C2. 

The derived ages for CI and C2 are also quoted in Table 1. In the next two sections we 
determine the luminosity and mass functions for both clusters and evaluate their slopes. 

4 Luminosity and Mass Functions 

4-1 Cluster Luminosity Functions 

In Figure 9 we show completeness corrected luminosity functions for CI and C2. Field 
contamination was eliminated by subtracting the field LF from that within the cluster region. 
The field LF was taken to be the average over several control regions equidistant from the 
cluster. The cluster regions used for determining their LFs were circles of 8" radius centered 
on each cluster. 

CI seems to have a slightly steeper LF than C2. LF slopes were obtained from linear fits 
to the points. We obtained 7 = Alog<^{Mv)/ AMy = 0.19 ± 0.03 for CI in the range My < 6 
and 7 = 0.12 ± 0.05 for C2 in the range My < 4. The best fit lines are shown in the figure. 
Upper limits were not included in the fits. Given the small number of LF bins and the low 
contrast of the clusters, the slope differences are not significant. In fact, the two LF slopes 
are similar in the common range used for the fits {My < 3). C2 LF basically ends beyond 
that while CI LF steepens. The 7 values are in agreement with that inferred from the work of 
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Flower et al. (1980) for NGC1868, a rich LMC cluster with similar age CI and C2 but larger 
mass. On the other hand, the inferred values for 7 are smaller than those typically fitted to 
younger globular clusters and to stellar associations in the LMC (Vallenari et al. 1993, Will et 
al. 1995a,b), 

In Figure 10, we show CI LFs for 3 radial bins; the two first are 1 pc wide, the last is 2 
pc wide. All 3 LFs shown are field subtracted and were scaled to the entire cluster area. The 
number of stars per unit area decreases with radius, as expected. The outermost ring contains 
only 3 bins in My with numbers significantly above the background. This ring includes the 
stars that make up the excess surface brightness seen beyond R = 2 pc in panel 6b. Despite 
the uncertainties, the LF looks shallower in panel 10c than in 10b, consistently with figure 6, 
where only an excess of light, not stars, is seen. In the inner rings, on the other hand, the LF 
becomes steeper with radius, providing evidence for mass segregation within CI. Thus, the 
bright stars beyond i? ~ 2 pc are likely to be background stars, belonging either to the general 
LMC field or to some stellar association superposed to CI. In fact, CI is situated in a rather 
messy border region between different clusters and associations, among them NGC2050, LH96, 
DEM261 and NGC157 (see catalogs by Lucke & Hodge 1970, Davies et al. 1976). The best fit 
slopes for the LFs in the two inner rings are 7 = 0.09 ± 0.04 and 7 = 0.25 ± 0.02. The outer 
ring lacks enough points for reliable fits to be made. 

It was not possible to split C2's LF into radial sectors, given its smaller contrast with the 
background. 

4-2 Cluster Mass Functions 

Comparing the MFs of clusters with different masses, ages, metallicities or environments 
can contribute to the debate about the universality of the initial mass function (IMF) and its 
evolution. 

The MF can be derived from the LF with the aid of a mass-luminosity (M-L) relation, 
which usually depends on metallicity and is often uncertain. However, the M-L relation is 
reasonably well known for the typical metallicity of the LMC, Z ~ 0.008 {[Fe/H] ^ —0.3), in 
the range of luminosities covered in this work. The mass functions for CI and C2 are shown 
in Figure 11. The M-L conversion was based on the Yale isochrone that best fits the clusters 
CMDs. Field subtraction proceeded in the same way as with the LFs. For field stars we used 
the same M-L relation as for the clusters. The derived mass of the evolved stars will obviously 
be in error, but these are subtracted off with the field, having little or no effect on the cluster 
MFs shown. 

Linear fits were made to the data points leading to F = Alog(f){M) / AlogM = — 1.8 ± 0.7 
and — 1.2 ± 0.4 for CI and C2, respectively. For C2, however, a single power-law fit to the MF 
is inappropriate, since its MF seems to be steeper for M 2 Mq. 

The derived slopes are subject to several sources of error. Isochrones with different 
metallicities or based on stellar models which incorporate convective core overshooting would 
change the M-L relation and therefore the MF slopes. The effect, however, is known to be 
small: ~ 0.1 in F (Elson et al. 1989, Sagar et al. 1991). Unresolved binaries lead to an observed 
MF which is flatter than the single star one. Sagar & Richtler (1991) investigated this issue 
and concluded that the amplitude of the effect is a function of the binary fraction and the 
MF slope itself. The observed slopes derived in this paper may be consistent with a true 
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r ~ — 2 and —1.5 for CI and C2, respectively, if the fraction of binaries is as large as 0.5. 
Finally, completeness errors may distort the actual shapes of LFs and MPs, especially if mass 
segregation is present, as seems to be the case in CI. Given the small size of the clusters, it was 
impossible to split them into many annuli for which completeness and mass functions could be 
measured independently. Given all these uncertainties we cannot rule out a single MF shape 
accounting for both of them. 

The LP and MF slopes are listed in Table 2, along with estimates of total luminosities 
(Mv), masses (M) and mass-to-light ratios (M/L). These latter were obtained by adding up 

the contributions of all non-saturated stars with My ^ 6.5 (M ^ 0.9 Mq). Thus, the absolute 
magnitude and mass estimates listed should be considered as lower limits. 

The resulting MF slopes are bracketed by most found in the literature. These were all 
derived for clusters much larger than the ones studied here. Elson et al. (1989), for instance, 

obtained — l.S'^F'^ — 0.8 for six rich young (r ~ 30 — 50 Myrs) LMC clusters in the mass 
range 1.5 - M/Mq - 6. Such flat IMFs were in large contrast with the results of Mateo 
(1988), who found slopes (—2.5 ^ F — 4.6) for six LMC clusters spanning a wider range of 
ages and metallicities. Sagar & Richtler (1991) derived intermediary slopes, F ^ —2.1, for five 
LMC clusters within the mass range 2 < M/Mq < 14 and spanning an age range of 10 — 100 
Myrs. 

More recently and using HST data. Hunter et al. (1995) obtained F = —1.22 for R136, 
a ^ 2 X 10^ Mq, very young and compact system in the center of 30 Doradus, in the range 
2.8 < M/Mq < 15. Our results suggest flat MF slopes for low mass LMC clusters as well. 

A more meaningful comparison would be with open clusters in the Galaxy. Ground- 
based studies include that of Sagar et al. (1986) who derived MF slopes for 11 open clusters. 
One of them (N1778C) has an age comparable to those of CI and C2 in §3.2. Its slope 
is F = —2.04 ± 0.44 in a mass range comparable to ours. This is consistent with CI and 
marginally consistent with C2. Other more recent ground-based works have led to shallower 
slopes, F — 1, but usually in the mass range M ^ 0.5 Mq (Comeron et al. 1996, Hambly et 
al. 1995 and references therein). 

4-3 Internal dynamics 

We now use the integrated properties and structural parameters listed in Tables 1 and 
2 to assess the importance of dynamical effects such as two body relaxation and tidal forces. 
Knowledge of the internal dynamical state of CI and C2 may help determining the extent to 
which such low mass clusters could contribute to the field population. 

Using the LMC tidal field derived by Elson, Fall & Freeman (1987), we estimate CI and 

C2 to have similar tidal radii, in the range 6 ^ Vt (pc) ^ 15, implying rt > fyis (see Table 1). 
Therefore, the clusters are so far not tidally truncated. Both, however, may have experienced 
stronger tidal fields or had close encounters with other more massive clusters or associations 
during their ^ 10^ years of existence. 

Is two-body relaxation relevant to the internal dynamics of CI and C2? The relaxation 
time scale for stars with a mass M at a distance r from the centre is given by 
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where v is the typical velocity of such a star within the cluster. Assuming hydrostatic 
equilibrium in the central regions, the central dispersion velocity in the cluster is given by 

a\0) r^Gpo rl (2) 

where Vc is the core radius and po is the central mass density. Using the relation between 
projected and spatial densities given by Djorgovski (1993), we estimate po ^ 16 Mq/pc^ for 
CI and Po ^ 6 Mq/pc^ for C2. These values for the central densities imply cr(0) ~ 0.4 km s"-*^ 
for CI and cr(0) ~ 0.2 km s"-*^ for C2. Inserting these velocities into equation (1) we infer that 

trei ^ 100 Myrs for a 1 Mq star in the central regions of both clusters. Therefore, both CI 
and C2 should be relaxed in their central parts. The evidence for mass segregation in CI is 
consistent with that. Dispersion velocities just about twice the values inferred would disrupt 
the clusters in ^ 50 Myrs. Given their estimated ages, CI and C2 are probably bound. 

5 DISCUSSION AND CONCLUSIONS 

We presented V and I photometry of two WFPC2 fields located near the eastern end of 
the LMC bar. Stars as faint as My ~ 7 (M ~ O.SMq) were detected in each of them. Each 
field contained one small open cluster, implying a large number density of such systems in the 
LMC. 

CI is roughly symmetrical in shape, has a steep density profile and a mass function slope 
of r = —1.8 ± 0.7. From isochrone fits, we infer an age of r ~ 200 — 500 Myrs for it. Its 
estimated mass is ~ 1500 Mq and its absolute magnitude. My ~ —3.5. The derived values for 
the luminosity, mass, MF slope and age are mutually consistent in light of the recent stellar 
population synthesis models of Girardi et ah (1995). These authors use the photometric models 
of single stellar populations calculated by Bertelli et al. (1994) in order to revise the relation 
between integrated photometric properties, age and metallicity of LMC clusters. From their 
Figure 13, we infer that a M ~ 2 x IO^Mq cluster with a MF slope a bit shallower than a 
Salpeter one (F = —2.35; Salpeter 1955) would have My ~ —3.5 at an age r ~ 10^ years. 

Marginal evidence of mass segregation within r ^ 2 pc from Cl's centre was found. This 
is consistent with the short relaxation time scale expected for M 1 Mq stars in its central 
regions. Tidal effects from the LMC should not be relevant; we derive r^i ^ 1 pc from Cl's 
surface brightness profile, which places its member stars well inside the estimated tidal radius 
(8 ^ Vt (pc) ^ 15). Given the central density and size estimates for CI, its stars should have a 
maximum central velocity dispersion of ~ 0.4 km s~^. It would be interesting to confirm that 
with observations. 

C2 seems less massive but just as luminous as CI (M ^ 4OOM0, My ^ —3.5). This 

implies fiatter luminosity and mass functions. We obtain F = —1.2 ±0.4 for C2. For My ^ 3.5, 
C2's LF is similar to Cl's but it drops off for fainter magnitudes. C2 is more irregular than 
CI and has shallower density and surface brightness profiles. We could only derive an upper 
limit of r 500 Myrs to the age of C2 based on its CMD due to saturation and reddening 
uncertainties. Based on the results of Girardi et al. (1995) and our mass. My and F estimates, 

we would obtain more stringent limits: r ^ lOOMyrs. C2's size was harder to quantify, given 
its fiat profile. However, just as CI, it is also likely to be contained within its tidal radius and 
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to have undergone significant core relaxation: tj-ei ~ 100 Myrs for M ^ IMq stars within the 
central 2 pc. 

The chance discovery of two small clusters in two MDS fields within the LMC suggests 
that these could be more common in the LMC bar or 30 Doradus region than previously 
anticipated. Such systems may have significantly contributed to the field star population if 
they were even more common in the past and got disrupted by larger clusters or by the LMC 
tidal field. Even though that does not seem to be the case for CI and C2, other similar systems 
may have been subjected to stronger tidal fields from both the LMC or its bar or from more 
massive clusters or associations. Alternatively many such clusters might have been unbound 
since their birth, their member stars streaming away after a few hundred Myrs. It would 
interesting to confirm or not the existence of a large population of low mass and luminosity 
clusters in the LMC with other observations using the high resolution of HST. 
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Table 1. Clusters structural parameters. 



Parameter 


CI 


C2 


(To (stars/pc^) 


30 ±8 


8±6 


Ra (pc) 


- 2 


- 1.5 


//y(0) (mag/arcsec^) 


20.2 ±0.4 


20.3 ±0.5 


(pc) 


- 2 


- 2 


^/iz (pc) 


- 0.9 


- 0.7 


Dyis X C?i,js (pc) 


2.9 X 2.4 


4.6 X 3.9 


T{Myrs) 


200-500 


- 500 
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Table 2. LF and MF slopes. 



Parameter 


CI 


C2 


7 = Alog^{Mv)/AlogMv 


0.19 ±0.03 


0.12 ±0.05 


T = Alog(j){M)/AlogM 


-1.8 ±0.7 


-1.2 ±0.4 


Ml 


3.5 


- -3.5 


Mass^ (Mq) 


- 1500 


- 400 


M/L (Mq/Lq) 


- 0.7 


- 0.2 



1 Due to stars in the range -0.5 ^ My ^ 6.5 (0.9 ^ M/Mq ^ 4.5) only. 
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FIGURE CAPTIONS 



1- The sharpness (s) and parameters, as obtained by fitting each detected object in chip 
4, field 1 to the WFC2 psf, are shown as a function of Feoe magnitudes, a) s. b) x^- 

2- a) CMD for all 3 WFC2 chips in field 1. The magnitudes and colours are in the HST 
WFPC2 system. The data are corrected for aperture and reddening effects. Superposed 
are Yale isochrones corresponding to Z=0.01 stars with ages 0, 200, 500, 1000 and 2000 
Myrs. b) The same as in a but now for the 3 WFC2 chips in field 2. 

3- a) Veoe completeness functions for field stars in field 1, chip 4 (squares and dashed line) 
and in field 2, chip 3 (crosses and dotted line). Error bars are standard deviations from 
different realizations of ADDSTAR, as explained in the text, b) Same as in a but now 
cluster completeness functions are shown: squares and dashed line: CI; crosses and dotted 
lines: inner and outer C2. Error bars are shown only for one of the cluster completeness 
functions to avoid confusion. 

4- WFPC-2 image of LMC field 1 cluster, CI. The whole field shown is about 50" on a side. 

5- WFPC-2 image of LMC field 2 cluster, C2. The whole field shown is about 40" on a side. 

6- a) Stellar density profile within the region occupied by CI. The upper curve gives the 
surface density of stars as a function of distance from the cluster centre, without any 
correction for field contamination. The points indicate the background corrected profile. 
Arrows correspond to upper limits. The adopted background density level is given by the 
solid horizontal line. The dotted horizontal lines show the background densities estimated 
separately for each WFC2 chip, b) Y band surface brightness profiles and background 
levels corresponding to the number density ones shown in panel a. 

7- a) The same as in panel 6a but now for C2. b) The same as in panel 6b but now for C2. 

8- a) The colour magnitude diagram of stars within CI region. Yale isochrones of 0, 200, 
500, 1000 and 2000 Myrs have been superposed to the data. The dashed line indicates the 
saturation level, b) The same as in panel a but now for C2. c) CMD for field comparison 
region close to CI. d) CMD for field comparison region close to C2. 

9- The luminosity functions of CI and C2 are shown as solid and open squares, respectively. 
The error bars incorporate Poisson fluctuations as well as uncertainties in the complete- 
ness corrections and in the background subtraction. Upper limits are shown with arrows. 
The lines show power-law fits to the data. 

10- The luminosity function for CI in different radial annuli, as indicated. The numbers 
in all panels are scaled to the total cluster area. Error bars again incorporate statistical 
uncertainties as well as completeness and backgroimd errors. Arrows correspond to upper 
limits. The solid line shown in the two upper panels is the best fit to a power-law. 

11- Mass functions of CI (solid symbols) and C2 (open symbols). The error bars incorporate 
Poisson fluctuations as well as uncertainties in the completeness corrections and in the 
background subtraction. The lines show the best flt power-laws in each case (CI: solid 
line; C2: dashed line). 
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